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ABSTRACT 

In this Letter, we present sensitive millimeter SiO (J=5-4; v=0) line observations of the outflow arising from 
the enigmatic object Orion Source I made with the Atacama Large Millimeter/Submillimeter Array (ALMA). 
The observations reveal that at scales of a few thousand AU, the outflow has a marked "butterfly" morphology 
along a northeast-southwest axis. However, contrary to what is found in the SiO and H2O maser observations 
at scales of tens of AU, the blueshifted radial velocities of the moving gas are found to the northwest, while the 
redshifted velocities are in the southeast. The ALMA observations are complemented with SiO (J=8-7; v=0) 
maps (with a similar spatial resolution) obtained with the Submillimeter Array (SMA). These observations also 
show a similar morphology and velocity structure in this outflow. We discuss some possibilities to explain these 
differences at small and large scales across the flow. 
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1. INTRODUCTION 

Since its discovery, the radio Source I in the middle 
of the Kleinmann-L ow nebula in Orion, has always been 
an en igmatic object (Menten & Reidl l 1 995t IChurchwell et afl 
1987). Recent proper motion measurements suggest that the 
Orion Source I might have been part o f a multiple system 
that disintegrated some 500 years ago dGomez et all 12008: 
Rodrig uez et al.ll200H iGoddi et alJl20Ul) originating an ex- 
plosive outflow in the Orion BN/KL region dZapata et alJ 
2009; [Bally et alJ l20Tlt iZapata et alJ 1201 lah . IBallv et alJ 
d201 11) proposed that probably the Orion Source I is a short- 
period binary formed just before the stellar ejection and the 
formation of the explosive outflow. It is estimated that the sep- 
aration of the bin ary is only a few Astronomical Units (AU) 
(IBallv etall2011h . 

Strong H2O and SiO vibrationally excited (v = 1 and 
2) masers are concentrated very close to the Orion Source 
I. The H2O masers are distributed in an elongated pattern 
along a northeast-so uthwest axis and seem to be expand- 
ing about Source I dGenzel et all 119811) . The bulk of the 
SiO masers are found surrounding Source I and form an X- 
shaped structure; the emission in the southeast arms is pre- 
dominantly blueshifted, while the emission in the northwest 
is predominantly redshifted. The intermediate-velocity emis- 
sion is observed connecting the red and blue sides of the 
emission distribution. The SiO masers seem to be arising 
from a wide-angle bipolar wind emanating from a rotating, 
ionized edge-on disk with a southeast-northwest orientation 
(Matt hews et al.|[2010l) . The ionized disk has been imaged at 
7 mm. and appears to be surround ing a young massive star 
dReid et al]l2007T: IGoddi et alJl20ll . 

At scales of a few hundreds of AU, iPlambeck et all (120091) 
imaged the outflow from the Orion Source I using CARMA 
(The Combined Array for Research in Millimeter-wave As- 
tronomy) in the J = 2 — 1 from the vibrational ground state 
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of SiO, and found that the central velocity channels showed 
the limb-brightened edges of two cones centered on Source I. 
However, since the outflow appears to lie nearly in the plane 
of the sky, there is no measurable red-blue velocity asymme- 
try in the cones, as observed for the masers. 

We report ~ 1.5" resolution images in the SiO J=5-4; v=0 
obtained with the Atacama Large Millimeter/Submillimeter 
Array 4 (ALMA) from the Orion Source I. These observations 
are complemented with SiO (J=8-7; v=0) maps obtained with 
the Submillimeter Array 5 (SMA). Both observations reveal 
that the outflow from the Orion Source I has a "butterfly" mor- 
phology along a northeast-southwest axis. The blueshifted ra- 
dial velocities are found to the northwest, while the redshifted 
velocities toward the southeast. 

2. OBSERVATIONS 

2.1. ALMA 

The observations were made with sixteen antennas of 
ALMA on January 2012, during the ALMA science verifi- 
cation data program. The array at that point only included 
antennas with diameters of 12 meters. The 120 independent 
baselines ranged in projected length from 18 to 253 kA. The 
phase reference center for the observations was at a 72000.0 = 
05 h 35 m 14. s 35, £72000.0 = -05°22'35'.'00. The primary beam of 
ALMA at 217 GHz has a FWHM ~ 30". The emission from 
the outflow powered by the Orion Source I falls very well in- 
side of the FWHM. 

The ALMA digital correlator was configured in 4 spectral 
windows of 1875 MHz and 3840 channels each. This pro- 
vides a spectral resolution of 0.488 MHz (~ 0.7 km s" 1 ) per 
channel. As a special case, this observation consisted of 5 
different spectral tunings resulting on 20 spectral windows. 
The spectral window where the line SiO (J=5-4; v=0) at a 

4 ALMA is a partnership of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA 
(Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Ob- 
servatory is operated by ESO, AUI/NRAO and NAOJ. 

5 The Submillimeter Array (SMA) is a joint project between the Smithso- 
nian Astrophysical Observatory and the Academia Sinica Institute of Astron- 
omy and Astrophysics, and is funded by the Smithsonian Institution and the 
Academia Sinica. 
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Fig. 1 . — SiO(5-4) (red line) and SiO(8-7) (blue line) spectra from the Orion BN/KL 
region obtained with the ALMA and SMA observations, respectively. The spectrum 
produced from the ALMA data is multiplied by a factor of two for convenience. In 
the ALMA data 1 Kelvin is equivalent to 9.04 Jy, while in the SMA data 1 Kelvin is 
equivalent to 1.73 Jy. 



rest frequency of 217.10498 GHz was detected, has a central 
frequency of 217.225 GHz. 

Observations of Callisto provided the absolute scale for the 
flux density calibration while observations of the the quasar 
J0607-085 (with a flux of 1.4 Jy) provide the gain phase cal- 
ibration. The total on-source integration time was about 100 
min. This time was divided to observe the 5 spectral tun- 
ings. For time-dependent gain calibration, the nearby quasar 
J0607-085 was observed every 7 minutes. 

The data were calibrated, imaged, and analyzed using the 
Common Astronomy Software Applications (CAS A). To an- 
alyze the data, we also used the KARMA software dGoochl 
1996). The resulting r.m.s. noise for the line images was about 
100 mJy beam -1 at an angular resolution of 1 ''78 x 1"23 
with a P. A. = -4.4°. We used uniform weighting to obtain 
a slightly better angular resolution. 

2.2. SMA 

The observations were obtained from the SMA, and were 
collected on December 2009, when the array was in its com- 
pact configuration. The 21 independent baselines in the com- 
pact configuration ranged in projected length from 16 to 70 
LI. The phase reference center for the observations was at 
a/2000.0 = 05 h 35 m 15M5, 5 /2 ooo.o = -05°22'35'.'0. Two fre- 
quency bands, centered at 335.3306310 GHz LSB (Lower 
Sideband) and 345.3306310 GHz USB (Upper Sideband) 
were observed simultaneously. The observations were made 
in mosaicing mode using half -power point spacing between 
field centers and thus covering the entire BN/KL region. The 
primary beam of the SMA at 345 GHz has a FWHM ~ 30". 

The SMA digital correlator was configured in 24 spectral 
windows ("chunks") of 104 MHz and 128 channels each. This 
provides a spectral resolution of 0.812 MHz (~ 0.7 km s ) 



per channel. The zenith opacity (T230Gffj)> measured with the 
NRAO tipping radiometer located at the nearby Caltech Sub- 
millimeter Observatory, was about 0.07 , indicating excellent 
weather conditions during the observations. Observations of 
Uranus provided the absolute scale for the flux density cali- 
bration. The gain calibrators were the quasars J0541-056 and 
J0607-085, while 3C454.3 was used for bandpass calibration. 
The uncertainty in the flux scale is estimated to be between 15 
and 20%, based on the SMA monitoring of quasars. Further 
technical descri ptions of t he SM A and its calibration schemes 
can be found in lHo et all (120041) . 

The data were calibrated using the IDL superset MIR, orig- 
inally developed for the Owens Valley Radio Observatory 
(OVRO. IScoville et al.|[l99l and adapted for the SMA. 6 The 
calibrated data were imaged and analyz ed in the sta ndard 
manner using the MIRIAD and KARMA ( iGoochll 19961) soft- 
wares. We set the ROBUST parameter of the task INVERT 
to to obtain an optimal compromise between resolution and 
sensitivity. The r.m.s. noise in each channel of the spectral 
line data was about 170 mJy beam -1 at an angular resolution 
of 2'.'9 x L'9 with a P.A. = +16.6°. The SiO J=8-7; v=0 was 
detected in the USB at a rest frequency of about 347.3 GHz. 

3. RESULTS 

In Figure 1, we show the total power SiO spectra of the 
ALMA and SMA observations toward the central part of the 
BN/KL region. The total area from the spectra was obtained 
covers approximately one arcmin 2 centered at the position of 
the Source I. No obvious maser emission is detected at these 
transitions. The spectra are very similar, but with that ob- 
tained at higher frequencies (SMA) being much stronger than 
the obtained with ALMA. This has also been observed in out- 
flows associated with young low mass stars (e.g. HH 211; 
iLee et al]l2007t iPalau et alJl2006h and is attributed to the SiO 
excitation within the outflow. The spectra have two central 
strong peaks, one peaking around of +6 kms 1 and the sec- 
ond approximately at +12 kms -1 . These peaks are associ- 
ated with the emission from the outflow of the Orion Source 
I around the systemic velocity of the cloud (+9 kms -1 ). The 
dip at 9 kms 1 could be self-absorption. The high-velocity 
structure has multiple components, some of them come from 
the outflow of the Orion Sour ce I, but others come from the 
molecular bubble reported by Zap ata et al.l (1201 lbb . and ob- 
jects located in BN/KL. We estimated that the radial velocity 
of the Orion Source I is about +6 kms' 1 , in g ood agree- 
ment with that value found by [Flambecket al. (200§|) and 
Matt hews et all d2010l) . 

The channel velocity map of the SiO (J=5-4; v=0) emission 
from the Orion Source I mapped with ALMA is shown in Fig- 
ure 2. The range of radial velocities (-15.5 t o +24.5 kms" 1 ) 
is in go od agreement with the one obtained by Plambeck et al. 
( 2009) shown in their Figure 1 . However, our maps with lower 
angular resolution seem to recover some extended parts of 
the outflow lost in their maps. This is more evident in our 
channels maps with central velocities of -15.5 kms~' and 
+ 14.5 kms -1 . In these channels, we observe some north- 
south structures with a "wing" shape, extending in an east- 
west orientation. The emission arising from the central veloc- 
ity c hannels could not be r esolved with our actual resolution, 
as in lPlambecket all (120091) . 

Integrated intensity maps (moment 0) of the SiO emission 

6 The MIR-IDL cookbook by C. Qi can be found at 
http://cfa-www.harvard.edu/~cqi/mircook.html 



05 3515.0 14.6 14.2 13.8 



RIGHT ASCENSION (J2000) 



Fig. 2. — ALMA channel map of the SiO (J-5-4; v= 0) line toward the Orion Source I. Every velocity channel is 5.0 km s 1 wide and is centered close to the position of the Orion 
Source I, Each box is labeled wi th the central LSR velocity in km s _I . A pink star marks the location of radio Orion Source I. The line represents the orientation of the outflow at a 
P.A.- 55° iPlambeck et aU2009l) . The contours are from -5, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, and 60 times 0.14 Jy beam" 1 , the rms noise of the image. The half-power contour 
of the synthesized beam is shown in the bottom left corner of the first panel. 
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Fig. 3. — Upper panels: ALMA (left) and SMA (right) integrated intensity (moment 0) color-scale and contour maps of the SiO thermal emission from the outflow arising from the 
Orion Source I, The blue color-scale and contours represent blueshifted gas with LSR velocities ranking from -15 km s -1 to —5 km s . The red color-scale and contours represent 
redshifted gas with LSR velocities ranging from +15 km s -1 to +25 km s -1 . The green color-scale and contours represent redshifted high velocity gas with LSR velocities ranging from 
+26 km s _1 to +30 km s _1 . In the left panel, the blue and red contours range from -10% to 90% of the peak emission, in steps of 5%. The emission peaks for the blueshifted and 
redshifted emission are 28 and 59 Jy beam -1 km s -1 . The green contours range from -25% to 90% of the peak emission, in steps of 10%. The emission peak for the redshifted emission 
(associated with the green color) is 27 Jy beam -1 km s . In the right panel, the blue and red contours range from -20% to 90% of the peak emission, in steps of 5%. The emission 
peaks for the blueshifted and redshifted emission are 140 and 350 Jy beam -1 km s -1 , respectively. The green contours range from 40% to 90% of the peak emission, in steps of 10%. 
The emission peak for the redshifted emission (associated with the green color) is 100 Jy beam -1 km s -1 . In both p anels, the synthesiz ed beams of the observations are shown in the 
lower right corner. Th e white dot marks the position of the Orion Source I, as derived from VLA proper-motion data (Gomez et al. 2008). Lower panel: This image was obtained from 
Matthews et al. 1 2010) and displays the velocity field of the 28 SiO v = 1 and v = 2, J = 1-0 emission surrounding the Orion Source I. 



obtained with ALMA and SMA are presented in Figure 3. In 
these maps, we only integrated over redshifted and blueshifted 
radial velocities associated with the outflow from the Orion 
Source I. The blueshifted gas emission was integrated over 
LSR velocities ranging from -15 km s _1 to -5 km s _1 . The 
redshifted gas emission was integrated over the LSR veloci- 
ties range from +15 km s _1 to +25 km s _1 . The very high 
redshifted velocities range from +26 km s _1 to +30 km s _1 . 
These maps reveal that the outflow from the Orion Source 
I has a marked "butterfly" morphology along a northeast- 
southwest axis, with the blueshifted radial velocities found to 
the northwest, while the redshifted velocities are toward the 
southeast. This velocity structure is differen t of what is ob- 
served at smaller scales, see the Figure 7 of Matthews et al. 
(2010) or our Figure 3. The ranges of radial velocities of the 
SiO masers are similar to those obtained by us. The wings 
of the "butterfly" shape extend several arcseconds along the 



east-west direction. 

At the middle of the "wings", we find very high velocity 
emission (see Figure 2, the green color-scale and contours) 
that also appears to emanate from the Source I, this compo- 
nent is observed at sm all scales with the same orientation by 
Matth ews et alj ( 120101) . In the ALMA map, with a better sen- 
sitivity and angular resolution than the SMA one, the tail ex- 
tends a few more arcseconds. The "butter fly" morphology can 
also be seen in the 0.5" resolution map of Wright et al. (1995) 
in the SiO (J =2-1; v=0) and in Figure 3 of Plambecket al. 
(2009) which was derived from 1' resolution data . The wings 
of the butterfly also show some differences between both 
maps. The red wing is only half, different from the blue wing 
that is shaped like an U in the ALMA map. Something seems 
to suppress the other red half. Maybe it is due to some density 
inhomogeneities in the surroundings of the bipolar outflow. 

4. DISCUSION 
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that the gas ejected ballistically at small scales from a rotating 
structure will, at large scales, appears as an expanding struc- 
ture, with the largest radial velocities overlapping at the center 
of the structure. This is not what is observed and we rule out 
this possibility. 




Fig. 4. — The central dot marks the position of the star and the inner circle indicates a 
rotating disk. Gas ejected ballistically from the disk will appear, at larger scales, radially 
expanding (and not rotating). 

In the following, we will discuss some possibilities trying 
to explain the nature of the velocity structure and morphology 
of the outflow emanating from Source I. 

4. 1 . Two bipolar outflows 

A likely possibility is to have two bipolar outflows ema- 
nating from a short-period binary star located in the middle 
of the ionized disk of the Source I. This is in agreement with 
what was proposed by Ba lly et al.l ([201 1 ) about the binarity of 
the Source I. One outflow could be moving in a nearly north- 
south orientation with its blueshifted side in the north and the 
redshifted side toward the south. The second bipolar outflow, 
would have an east-west axis, with its blueshifted side to the 
west, while its redshifted side is to the east. The difference 
in size of the components could be attributed to the molecular 
cloud inhomogeneities. The fact that a binary could have a 
strong precession moving in short time scales, might explain 
the difference in the redshifted and blueshifted velocities ob- 
served at different scales. The SiO masers could be tracing 
the innermos t parts of the two bipo lar outflows ejected at the 
present time (Matth ews et al.ll2010t) . However, this hypothe- 
sis could not explain satisfactorily well the high velocity tail 
(with a southwest orientation) observed in the middle of the 
outflows, and the observations o f the single "cone" outflow 
found by Plam beck et al.l d2009i) at their central channel ve- 
locities. 

4.2. Ballistic ejecta from a rotating disk 

We also consider the possibility that the thermal SiO traces 
gas that was ejected ballistically in the past from the small- 
scale disk associated with Source I. However, Figure 4 shows 



4.3. A single expanding wide-angle NE-SW outflow 

Another possibility is having a wide-angle expanding out- 
flow that is emanating from the ionized, near ly e dge-on, ro- 
tating disk reported by Matthe ws et al.l {2010) and iReid et al.l 
(2007). In this case, the SiO masers would trace the inner- 
most parts of the outflow and the thermal SiO emission would 
be tracing the outer parts of the outflow. The difference in 
the orientation of the radial velocities at different scales could 
be attributed to that maybe in the past the rotation of the disk 
was in a different sense. This could have happened if a large 
amount of matter fell directly to the disk and changed the ro- 
tation sense. This third idea might explain most of the char- 
acteristics observed in Fi gure 3 and those observed at small 
scales by lMatthews et alJ d2010l) and lPlambeck etalJ (120091) . 
However, the idea of a disk inverting its rotation direction re- 
quires of unusual events during the evolution of the disk and 
deserves further theoretical and observational studies. 



5. CONCLUSIONS 

In this paper, we have reported SiO (J=5-4; v=0) and 
SiO (J=8-7; v=0) observations obtained with ALMA and the 
SMA. We found a marked "butterfly" morphology along a 
northeast-southwest axis in the outflow from the Source I. 
However, the blueshifted radial velocities are found to the 
northwest, while the redshifted velocities toward the south- 
east, which is contrary to what is observed at very small 
scales. We discussed some possibilities to explain this phe- 
nomenon, however, any of these hypotheses requires further 
sensitive, and high spatial resolution observations. 
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